, 


Y | gesE Me nam ES] | F ; - e - 3 É ã 
ae A “en 
t+ a o | 
do. instintos 
a. - RAS pa SE 
a ] E Pe o E .” 
68 bride mm e 
pa 


. 
0... 
teto! 
.. 


0%6840) 
ms o 
le teto resta tat ata R E D UTO R E S DA 


VARÕES DE AÇO DE ALTA RESISTÊNCIA 
PARA BETÃO ARMADO TIPO A. 40 den 
(«Heliaço» e «SNT 40») mei DE VELOCIDADE “ii: 
TENSÕES DE SEGURANÇA mus TIPO «PONTA DE VEIO» neto 
pp 5] aponta” mais de cem anos dn 
: rogresso Ptereteieitata 
Ra — 4000 kg/cm? (SISMOS) pd gd, Pa iii 
fabricantes «==: 
Literatura e Tabelas de Cálculo ingleses dis 
fornecidas gratuitamente it 
R. DE S. JULIÃO, 41, 4.º D. — LISBOA de AUTO-LUSITAN IA Ennis 
ECO Ã LFRE DO D UAR T E b - etnia 
Tel.: 36 7817-3237 06-369477 eetedo AV. DA LIBERDADE 73-77 LISBOA nt 


aum Bras! | CIMENTO SEC | 


coberturas 


perfis especiais 
vigas pré- DO 
-esforçadas ES » 
pavimento trief N ( Iê 


protecção de 


segurança 
Fábrica: 


materiais 
Outão -—- Setúbal 
novobra 


gg 


um 


Telef. 228824/5 


Escritório: 
RR, do Comércio, 5S6-3.º 


LISBOA 2 
Telef. 328201/2/3 


E 
pavimentos 
em betão pré-esforçado 

homologados pelo Lab. Nac. de Engenharia Civil) 


SOCIEDADE DE PRÉ-FABRICAÇÃO E OBRAS GERAIS NOVOBRA, LDA. 
AY. EST. UNIDOS DA AMÉRICA, 100 -5º D.º- TEL. 774832 / TTIGO3 - LISBOA 5 


NÚMERO 381 


Director 


Victor M. M. Pinto 


Administrador 


Victor Coias e Silva 


Corpo Redactorial 

João M. Pereira 

João Soares Alves 

Rui Coutinho Torrinha 


Silvério de Sousa Mendes 


Editor e Proprietário 


Associação dos Estudantes 


do Instituto Superior Técnico 


Redacção e Administração 


Avenida Rovisco Pais 
Telefone 77 5449 


Lisboa 1 


Assinaturas 
5 N.os 10 N.os 
70800 130800 


Continente e llhas  — 
Ultramar, Brasil 
80800 1405800 


90800 150300 


e Espanha 


| 


Estrangeiro 


Número avulso — 20$00 


Não se publica em Agosto 
e Setembro. 

Os artigos assinados são da 
exclusiva responsabilidade dos 


seus autores. 


Composto e impresso na Tipografia 
forge Fernandes, Lda. — Rua da 
Cruz dos Foiais 1058 — Lisboa 


TECNICA N.º 981 


OUTUBRO DE 196% 


SUMÁRIO 


41 


12-50 


CAPA, 


ANO XLIII 


a — 


VOLUME XNX 


PUBLICAÇÃO MENSAL 


LUIS MALTEZ — Trajectória dum sistema. 
Trajectory of a Sistem. 


A. M, Ro MARRECAS — Da teoria da fonte móvel 
num meio infinito a duas dimensões. 


About the theory or a moving heat source in a two dimen- 
sional infinite media 


LOBO FIALHO e MACHADO GONÇALVES — O efeito das 
variações térmicas sinusoidais nas estruturas 
laminares de comportamento elástico linear. 


Efect of sinusoidal thermic variations in concrete shells of 
linear elastical behaviour. 


VASCONCELOS PAIVA — O dimensionamento e a 
optimização de pavimentos pré-fabricados. 


Dimensionation and optimization of pre-fabricated floors. 


MADEIRA COSTA — Dimensionamento de secções 
de betão armado com compressão, em relação à 
rotura. 


Dimensionation of concrete sections submited to compression, 
in relation to rupture. 


ANTONIO VASCONCELOS — Algumas considerações 
sobre a forma dos acessórios de tubagem — Nota 
Técnica. 


Notas Informativas 


Elementos sobre a produção e o consumo de energia na rede 
eléctrica nacional. 


Bibliografia 


[ractor a vapor. Fotografia cedida pela Caterpillar Overseas S. A, 


ÍNDICE ALFABETICO DOS 


Pág. 

A. Johnson & Co nigia Ltd. XXVUHI 
MES: we 5 : à XV 
Asea Eléctrica . v, VII, XXIX e XXXI 
Auto-Lusitânia - 2.º dacapa 
Bethencourt Bros XXXI 
Cerâmica de Vale de Lobos, Lda. . XX 
Cimento Secil. . . «. . «o» 2.ºdacapa 
Companhia de Cimentos Brancos . VI 
Construções Técnicas . XVII 
PM Eras g E RS du Na XX 
Electro-Arco, Ltd.* . +. . ... i 5 H 
Empresa de Cinitiitós de Leitdá XXXIV 
XXXI 


Empresa Electro Cerâmica . 
Empresa de Produtos Asfálticos .. X 
Empresa de Sondagens e Fundações 

Teixeira Duarte, Lda. - 3.º da capa 
Empresa Técnica de Equipamentes 

Eléctricos . 
Empresa Ultramarina de Sondagens 


e Fundações . XVI 
Fanafel . .. ' EO ss XVI 
Fundações Frank i, tdi sea é É I 
Garagem Vitória . . . à & XXXI 
Guedes & Almeida, dá VII 
Harker Sumner . +... XXI 
Hola. san gas, amais ce 2.º da capa 


ANUNCIANTES 


Pág. 
INA XXVI 
Instron +... cc. XXX 
|O. T, . ES s XXVI 
OMS «uvcmszaz uma a yo XX 
J.J. Silva & Irmão . cc... V 
). F. Azevedo e Silva «2. ..c. IV 
j:. Roma + « Dn ESA 4 VIH 
Jayme da Costa, Ldi | XIX 
Johanna Keller :vsusmaas aci V 
LUSOTÉCNA XIV 
MAGUE XI 
Novobra, Lda. . 2.' da capa 
Química Hoechst XXIX 
ROO sema nmms sa XXVII 
Siemens XXI 
Soc. Const. Amadeu Gaudêncio XXVI 
Soc. Brown Boveri . . a 4.º da capa 
Soc. Michaelis de Vasconcelos . .. XH 
Sociedade Portuguesa Cavan . . ... VI 
Somaguê » cm vmoso XXXIV 
Sondagens e Ends A. Cavaco . 0! 
Sondagens Ródio, Lda. . . . ... XXVI 
Sopecate . . . .. E a H 
Standard Eléctrica . . . . .... IX 
DE ectumão kras ES XXHI 
do corsa: EA VII 
Wild Portugal, Lda. . . ..... 1 


TECNICA Nº 381 


SORUEAE SONSSÁNIES 


JEM a puD as 98 e 


Joenape gRRÊS TSE: 


ne 1 E E 1 2.4 E 
Sa 


MD oga ano nesios 


= = 


. O rh me = 


bs am ma aa” É ui a puitE 
- z ' ” « 


COMPANHIA CIMENTO TEJO duas baterias de silos e ensacagem - ALHANDRA 


Assente sobre Estacas FRANKI 


FUNDAÇÕES FRANKIL?A 


RUA BRAAMCAMP, 13-3.º — LISBOA 1 — TEL. P, P. C. 5341 12/3/4 — ENDEREÇO TELEGRÁFICO: FRANKI — LISBOA 
PRÉDIO SANTOS GIL, 6.º ANDAR -—- SALA 10 — TEL, 91167 —- CAIXA POSTAL 1397 


ENDEREÇO TELEGRÁFICO, FRANKI-LOURENÇO MARQUES 


TECNICA ! 


== 
” SOPECATE 


SONDAGENS 


' 


agi 


NMONHOZC 


A 


Ei 
N 


ESCRITÓRIO: RUA DO ARSENAL 146 


e ESTUDOS, PESQUISAS E CAPTAÇÕES DE AGUA 


e SONDAGENS 


* REBAIXAMENTO DE NIVEIS AQUÍFEROS 


| 


sopecate 


sondagens - fundações 


2º TELEFONES 34010-320208 LISBOA 


| | ELECTRO-ARCO, L.' 


ELÉCTRODOS + POSTOS + ACESSÓRIOS 
PARA A SOLDADURA ELÉCTRICA 


* 


MATERIAL APROVADO PELOS 
LLOYD'S REGISTER OF SHIPPING 


BUREAU VERITAS 
GERMANISCHER LLOYD 


* FUNDAÇÕES 
LISBOA 
Rua Silva Carvalho, 239 
Tels. 6836 49-68 48093 
FÁBRICA E LABORATÓRIOS 
Venda Nova — Amadora 
Tels. 9706 13-971733 
SONDAGENS E FUNDAÇÕES A. CAVACO, LDA, PORTO 
no 
R. RODRIGO DA FONSECA. 62, RC ESQ — LISBOA -1 Rua do Bolhão, 216 


Tel. 21277 - 394473 


TELEF. 56 1171/72/73) 


TECNICA II 


O novo Nível automático 
de engenheiro 


WILD NAK 2 


com circulo horizontal 


ja fimo, 


HEERBRUCGG 


Para nivelamentos de média e grande precisão, 
com micrómetro de lâmina de faces plano- 
-paralelas para nivelamentos de grande precisão. 


Erro médio para 1 Km de 
nivelamento duplo: + 1,5 mm 


Erro médio para 1 Km de 
nivelamento duplo com micró- 
metro: + 0,4 mm 


Luneta de 30 > de aumento 
com imagem direita de obser- 
vação 


Solicitem catálogo ou demonstração 
aos representantes exclusivos 


WILD PORTUGAL LDA. at ss + 


TECNICA IH 


un 
| qui 


CARACTERÍSTICAS 

modelos de 4 a 13 elementos 

ganhos de 5,8,10,11 dB 

captação de qualquer canal de 5a 
transformadores de impedâncias de 300/75 “d 


extensa gama de acessórios: cabo blindado coaxial, adaptadores de impedâncias, 
misturadores, tomadas, amplificadores de sinal, divisores de baixada etc. 


J.F DE AZEVEDO E SILVA & C:, Lº” 
Rua S. Francisco de Sales, 2 andar, Tel. 654165 — Lisboa 


TECNICA IV 


ns” 
ne a 


ES VV 


“aa o 
— Cd 4 nes .— 
o do 


pasa 550 esc. 


CARACTERÍSTICAS 1 Tubo de aço, articulado em qualquer posição, várias Incli- 
nações, [hcilmente trensportável, não ocupa espeço quendo 
fechado, prencheis de 1,20 X 0,80 de madeirs, em tola. 
Pintura de estufa. 


Em exposição, outros modelos de madeira e metálicos e todo o mete- 
rlal de consumo e equipamento de desenho, escritório e pepeleris 


Distribuidores exclusivos, em Lisboa, do estirador da Fábrica Olaio 


J. J. SILVA & IRMÃO, LDA. 


rua silvs cervalho 232-b » lisbos - telefone 652750 | 680671 


MARCA -A.E.1 


FUNDAÇÕES 


ESPECIAIS 


JOHANN KELLER 
TIS E 


TECNICA 


O EMPREGO DO CIMENTO BRANCO... 
permite acabamentos 
mais perfeitos, mais 


duradouros e muito 
mais económicos 


Estude a vantagem do emprego do 


CIMENTO BRANCO LUSO 


IBRA 


Consulte os distribuidores gerais 


No Sul: SCIAL —T. do Corpo Santo, 15 — Telef. 66 61 86 — Lisboa 
No Norte: SCIAL — R. do Bonjardim, 205 — Telef. 25779 — Porto 


POSTES DE BETÃO PARA 
LINHAS ELÉCTRICAS 
MOSAICOS E AZULEJOS 
COLUNAS PARA 
ILUMINAÇÃO PÚBLICA 
MANILHAS PARA ESGOTOS 


TECNICA MI 


WAaBCO 
NARA 


WESTIN'GHOUSE MOTO-SGRAPERS AUTO-GARREGADORES 


AIR BRAKE COMPANY 


t onstructon F qupment Division 


Lecria (Ulioos | Oido! 


CAPACIDADE DE 11, 21 E 31 JARDAS CÚBICAS 


NÃO NECESSITA DE TRACTOR PARA EMPURRE! 
UMA SÓ MÁQUINA... UM SO OPERADOR... 


Escarifica... Carrega... Transporta... Descarrega... Em terrenos de toda a natureza ! 


A.“ o 


) 


“ - » A f à 2 asd 4 
o . ã SM FATE 
Ava A ? * + AA E DE 


ee 


asi dia sá Eae EV ao ope pe“ Deo eia SR PA Doe DEDO pie DD DP E es 
Ema, * a Dr de per ter nd Aa à Gens ST 
e = ad RAE SETA PELA PR RES ass Es TA . Fiapo Pp ci iu és: AS E A Era ro gs: 

E E cito 4 = Lp. E o 4 3 “a = x É e 
SA a mp "2. PE SA FA TA” SPA sa mm ão " = - - a é dot * q 3 e 


DISTRIBUIDOR EXCLUSIVO 


GUEDES... ALMETAS LDA. 


Rua Áures, 181-2.º a LISBOA-2 — TELEFS. 3250 80 e 327845 
Através das suas organizações de: LISBOA + LUANDA LOBITO +» NOVA LISBOA e Agências de: PORTO, GUINE, 


S. TOMÉ E PRÍNCIPE e CABO VERDE fornecerá todos os elementos técnicos, para criteriosa escolha, assistência técnica e sobresselentes. 
Ara 


TECNICA VII 


EE RR SN 
e SOLDADURA 
ELÉCTRICA 


e MATERIAL 
DE PROTECÇÃO 


MARCA - A.E.2 


Aparelhagem 
de medidas 
eléctricas 


+ indicadora 
+ reguladora 


+ registadora 


J. ROMA, L.?* 


P.da Figueira, 12,1. P. do Município, 309-6.ºS 3 | 


LISBOA Tel. 865151 PORTO Tel. 0236732 


TECNICA VII 


na da 


as ul 


e day dE ay SS 
, 


material mais usado no isolamento térmico 
das construções 


vd SETH, LDA, .- Rua Filipe Folque, 10-1.º + Lisboa - Tel. 550156 


Tree meo + me 


PARA COLAGEM DE TACOS DE MADEIRA USE: 


Sh, DE Php 
Qf U, 


S 


Também proprio para: 
impermeabilizar terrenos, 


paredes, fundações, etc. FÁBRICA EM SA CA VÉM 


RUA FILIPE FOLQUE-10 ""— LISHA 
TELEF. 53 01 56 —TELEG EPA! JA-LISI' A 


IFOCNICA X 


O CRESCIMENTO COMO FACTOR DE DESENVOLVIMENTO 


Tal como as árvores, quando bem 
cuidadas, criam raizes, crescem, 
ramificam-se e dão bons frutos, 
assim as empresas, quando bem 
dirigidas, definem uma posição, 
desenvolvem-se, diversificam as 
suas actividades e lançam nos 
mercados produtos de alta qua- 
lidade. 

A Standard Eléctrica, pioneira da 
indústria das Telecomunicações 


e Electrónica em Portugal, que 
em 1932 não dispunha de mais 
de uma escassa dezena de em- 
pregados, conta hoje com mais 
de 2200 pessoas nas suas fábri- 
cas de Lisboa, Cascais e Luanda. 
O ritmo de crescimento e ex- 
pansão tem sido de tal modo 
elevado que o volume de ven- 
das passou de 30000 contos em 
1957 para 300 000 em 1967, pre- 


Standard Elecirica 


vendo-se que em 1970 a produ- 
ção total da Standard Eléctrica 
atinja o valor de 1500000 con- 
tos, dos quais 1000000 se des- 
tinam à exportação. 

A Standard Eléctrica é, pois, um 
valor positivo na indústria nacio- 
nal, que tem como base o ele- 
vado nível técnico e profissional 
e ainda o magnifico esforço de 
todos os que nela trabalham. 


TECNICA IX 


Comunicações 


equipamentos 


e sistemas completos 
* Cuidadosamente 
projectados 


« Amplamente 
experimentados 


* Completamente 
integrados 


Registador Sanborn 


HEWLETT PACKARD 


Registador digital 


aparelhagem electrónica Dymec 
de medida e registo 


Osciioscopios, Oscilógrafos, Voltimetros 
electrónicos, Geradores de sinal, 
Contadores electrónicos, Registadores 
digitais, Fontes de alimentação 


Oscilocópio HP 


Contador Electrónico MP 


VJELEGT ERA MESA TENCA DE COUPAMENTOS FLECTRICOS, SARL 


TECNICA XI 


Os para- 
raios 
Oerlikon... 


= = estão na van- 
guarda do 
progresso 
e protegem 
instalações 
em todo 
o mundo 


Descarga de longas linhas 
em vazio, protecção de 
postos completos, transforma- 
dores, alternadores ou 
consumidores ligados a linhas 
aéreas. Para cada função, 

o pára-raios apropriado. 


Corrente de 
descarga 5 kA 


380 & 500 V 


OERLIKON 


Ca a 


Corrente de descarga Corrente de descarga Insuflação magnética Para indicações mais 

5 e 10KkA 10 kA Ondas Corrente de descarga amplas consultar 
respectivamente rectangulares 10 kA, Ondas 

Ondas rectangulares 150A 2000 us rectangulares 98 nOSSus representantes 
75A 1000 us e Membrana 1000A 2000 us Sociedade Michaélis 

150A 2000 us anti-explosão Membrana anti-explosão de Vasconcellos, Ltda. 


Porto e Lisboa 


Ateliers de Construction Oerlikon - Zurich (Suiça) 


Representantes para Portugal é Ultramar: Sociedade Michaélis de Vasconcellos, Lda., POTIO E Lishoa 
TECNICA XII 


NUMERO 351 OUTUBRO DE 1965 


TRAJECTÓRIA 


RESUMO 


Apresenta-se uma possível forma de descrição matemá- 
tica da trajectória dum sistema, Este problema particular 
integra-se no breve estudo dos sistemas iniciado pelo autor 
na “Tecnica” n.º 378, com o artigo “Introdução ao Estudo 


>, 


dos Sistemas”, no qual devem ser procuradas as definições 


das noções básicas agora utilizadas. 


1. INTRODUÇÃO 


No nosso artigo “Introdução ao Estudo dos 
Sistemas”, publicado na “Técnica” n.º 378, 
chamámos, de forma simplista, trajectória dum 
sistema a qualquer sucessão de imagens num 
espaço cartesiano de observação. 

Sempre que não se registassem ciclos evolu- 
tivos, essas sucessões poderiam ser simples 
conjuntos totalmente ordenados de imagens, em 
que a relação de ordem definiria a anterioridade 
na sucessão. Se, porém, houver ciclos, em que se 
repitam imagens anteriormente verificadas, sere- 
mos forçados a imaginar outra forma de definir 
a anterioridade na sucessão. 

Assim, para dar generalidade à teoria, teremos 
de nos apoiar num conjunto totalmente ordenado, 
estranho ao sistema, que servirá de base de refe- 
rência e que, por correspondência com o espaço 
das imagens, permitirá destrinçar a ordenação 
destas na sucessão. 

É evidente, porém, que há certa arbitrariedade 
na escolha da base de referência, pelo que iremos 
encontrar diferentes tipos de correspondência 
conforme a base escolhida. 
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ANO NLIII 
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DUM SISTEMA 


por LUÍS MALTEZ 


Assistente do LS.T. 
SYNOPSIS 


A mathematical description of the trajectory of a aystem 
is presented, This is a particular problem soithin the 
brief study of systems started by the author in “Técnica” 
n.º 378, with the paper “Introdução av Estudo dos Siste- 
mas”. In that paper can be found the definitions of the 
basical notions here used, 


Por este motivo, para comparação do compor- 
tamento de dois sistemas, teremos sempre de 
referir ambos à mesma base. Do mesmo modo, 
se estivermos a estudar um único sistema, mas 
quisermos que esse estudo sirva para futuras 
comparações, teremos de especificar bem a base 
de referência, em condições de podermos recons- 
tituí-la para outros sistemas. 

Pelo contrário, a trajectória dum sistema, enca- 
rada como simples sucessão dos seus estados, 
ou das imagens destes, pode abstrair das peculia- 
ridades da base de referência. 

Dentro desta óptica, começaremos por estabe- 
lecer sequências de estados, em correspondência 
com uma base totalmente ordenada, mas poste- 
riormente definiremos uma relação de equivalên- 
cia entre sequências, de modo a permitir definir 
as trajectórias como classes de equivalência. 


Notas à simbologia 


Seja uma aplicação f de A em B. Seja 


FCAXB (1) 


o gráfico da aplicação f. Então f é a tríade 
f= (FAB) (2) 


O gráfico F define também uma relação 
binária 4 em A X B. 

Chamaremos f (a) a imagem do elemento 
a: A pela aplicação f, a qual é sempre um e um 
só elemento b: B. Usaremos o símbolo F (a) 
para representar o corte do gráfico F pelo ele- 
mento a: À, o qual é, neste caso, o conjunto 
(bj « P (B). Anâlogamente se representa pelo 
símbolo F (b) o corte de F por b:B, o qual é 
um sub-conjunto de A, podendo o seu cardinal 
ser qualquer. 

Para a imagem dum sub-conjunto S : P (A) pela 
aplicação f e pelo gráfico F, usaremos respectiva- 
mente os símbolos £f<S>eF<S>. 

Note-se que é sempre 


f<S>=E<S> (3) 
Espaço de observação 


Como se pode verificar no nosso artigo acima 
referido, se um sistema S for observado por um 
conjunto 


C=(Cjicl) (4) 


de critérios de observação, o espaço de observa- 
ção é o conjunto X V,, produto cartesiano dos 
jel 


conjuntos imagem por cada critério, eventual- 
mente organizado como espaço vectorial. 


A imagem do sistema S em X V éo conjunto 
j:l 


C<S>=([E:X V:(de:5) [E=C, (e)]) 


cl 


(5) 
Chamando 
tem-se ôbviamente 
5, E X V, (7) 


2. SEQUÊNCIAS 
Escolha da base de referência 


O problema que se põe no estudo dum sistema 
é o da previsão da sua evolução a partir dum 


certo estado conhecido. Por isso, as sequências 
de estados a considerar serão sempre sequências 
limitadas à esquerda. Quer isto dizer que pode- 
mos sempre usar como base de referência um 
conjunto totalmente ordenado com primeiro 
elemento. Resta analisar qual deve ser a sua 
potência. 

Se o sistema em estudo tiver um número 
finito ou numerável de estados, isto é, se 5, for 
finito ou mnumerável — o que não impede que o 
espaço cartesiano de observação x V, possa 

1 E 
ter potência superior — e se o sistema S evoluir 
sempre para um estado de estabilidade absoluta, 
então as sequências de estados serão sempre 
finitas, pelo que a base ordenada pode ser finita. 

Se, porém, sendo ainda S, finito ou numerável, 
o sistema contiver situações instáveis ou astá- 
ticas, isto é, se o sistema puder ter evolução 
divergente ou cíclica, já só poderemos descrever 
completamente as sequências sobre uma base 
numerável. 

Em todos os outros casos encontrados na prá- 
tica, verifica-se que ê sempre suficiente que a 
base de referência tenha a potência do contínuo 
sem último elemento, isto é, o tipo de ordem da 
semi-recta real não negativa. 

Como nada há de proibitivo na utilização duma 
base de potência superior à mínima necessária, 
vamos apenas tratar o caso geral, escolhendo 
para base de referência o conjunto totalmente 
ordenado (R +, <<) dos números reais não nega- 
tivos. Todos os casos particulares apontados se 
poderão obter por posteriores restrições das 
sequências obtidas em R * a um seu sub-conjunto 
numerável ou finito, ordenado pela ordem indu- 
zida. Veremos mais adiante como efectuar tais 
restrições. 


Sequências de estados 


Podemos obter uma sequência de estados asso- 
ciando a cada elemento da base de referência um 
e um só estado de S, ou melhor, um e um só 
elemento E : 5. 

Nestas condições, chamaremos sequência de estados 
a qualquer aplicação c de R* em S,. 

Ora, por definição de potência de conjuntos, o 
conjunto S/"* é o conjunto de todas as aplica- 
ções de R? em S,, pelo que podemos chamar 


R+ ; = 
a Si "o conjunto das sequências. 
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Cada sequência 7 : Sj8 + será então uma tríade 


c=(H,R*,S) (8) 
em que 


HCR'xs£S, (9) 


é o gráfico de 7. Cada elemento do gráfico H é 
um par ordenado do tipo (i, E) em que i ER" 
dá a localização na base de referência da fase 


evolutiva de S em que a imagem do seu estado 
em X V éEs 9, | 
jel 


3. SEQUÊNCIAS REDUZIDAS 


Partição de R' em intervalos conexos 


A “ R - “ 
Cada sequência 7:S|' permite operar uma 
partição da base R* em intervalos conexos, em 
cada um dos quais se mantém a imagem por 5. 


. “ 1: é je 
Para isso definiremos a família Oooes 8 de 


relações de equivalência em R”, como segue: 


E seca iê RA. as 
Dados i, je: R' e dado 7 : S/“”, diz-se que 
ijoj see só se 


(vk [(i=k<j)v(j<k<i) > (c(k)=c(i))] (10) 


Demonstra-se (Ver Anexo) que 7; é uma rela- 
ção de equivalência em R”, visto que goza das 
propriedades reflexiva, simétrica e transitiva. 

Podemos então determinar o conjunto cociente 
R*/ Ja de Rº pela relação %, que é o conjunto 
das classes de equivalência definidas em R* 
por 5. Cada uma destas classes (Ver Anexo) é 
um intervalo conexo em R', cujos pontos têm 
todos a mesma imagem por 7. 

É óbvio que podem resultar classes com um 
único ponto i « R*. Chamar-lhe-emos ainda in- 
tervalo — será o intervalo fechado [i, il. 


Base de referência reduzida 


O conjunto cociente R* / s acabado de formar 


é susceptível de ser totalmente ordenado por 
uma relação de ordem que simbolizaremos < e que 


se infere da relação <", como segue: 


6 
Dados 1, ,1I, R*/%,, diz-se que Ma I, 
se e só se 


(1, E L;) A (a, ' 13) Ê ha X Ig) (i, “ 14)] (11) 


| 
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Sendo 1» e I3 intervalos conexos em R* e mú- 
tuamente disjuntos, por serem diferentes classes 
de R*/%,, conclui-se que o quantificador exis- 
tencial q pode ser substituído, na definição da 
relação < , pelo quantificador universal V. 

Obtivemos deste modo um conjunto total- 
mente ordenado (R” /9s, <) de potência igual 
ou inferior à de (R*, <<); chamar-lhe-emos base 
de referência reduzida segundo 5. 

Repare-se que R” / Ja tem como primeiro ele- 
mento a classe que inclui o primeiro elemento 
de R” (o ponto zero). 


Sequência de estados reduzida 


A definição de R*/% permite-nos definir 
imediatamente, a partir de 7, uma nova aplicação 
:*, agora de R */Js em Sr, dizendo que 


do 
R a ad 
Dado qualquer 7 «: S/ |, chama-se sequência 
reduzida de 3 a aplicação 


ct = (H*, R*/%, S1) (12) 
de gráfico 
Ht=[(1,E)eRtSXS: [E] = 
= <1>) (13) 
Temos assim que 
(vie R*/%) (3<1I>= | (D|) (14) 


porque, sendo I uma classe de equivalência em 


R” segundo “Ys, todos os seus elementos k têm 
a mesma imagem q (k) e portanto o próprio I 
tem por imagem 7 < | > um conjunto com um 
único elemento q (k), isto é 


e<I>= |c(k)| (15) 


1 
. = R “ 
Em resumo, cada aplicação 7:S,  deuorigem 


à geração de uma base de referência reduzida 
G 


(R'/%5, =), sobre a qual foi definida uma apli- 
cação c* com imagens em S,. A imagem de 
uma classe por q* é a imagem de cada elemento 
dessa classe por 9. 


Seja Lp: o conjunto de todas as aplicações 
3 * assim obtidas. A forma utilizada na definição 


peca RT 
de 7* corresponde a uma aplicação v de Si 
sobre Zp+, a qual é bijectiva, como é evidente. 


Será 
0* = b (6) (16) 


o = q (9º) (17) 


. Eid R+ 
É assim equivalente falar na sequência 5:5, 
ou na sequência reduzida q: 2,1. 


4. TRAJECTÓRIAS 


A sequência reduzida de uma sequência 7 con- 
serva toda a informação sobre a correspondência 
inicialmente estabelecida com a base de referên- 
cia Rt. Com efeito, dado 


ot = (H*, Rt/%, 51) (18) 


conhecemos R*/7, e a imagem q* (I) de cada 
[e Rt/%,. Então, a sequência original q será tal 
que 


(yk =) (3 (k) =0* (1)) (19) 


É assim possível reconstituir 7 a partir de 7*, 
o que confirma que, como dissemos em 3,, a 
aplicação y é bijectiva. 

No entanto, a criação das sequências reduzi- 
das prepara o caminho para a abstracção da base 
de referência e, além disso, informa-nos já sobre 
o cardinal mínimo que deve ter essa base. Ele 


terá, com efeito, que ser pelo menos igual ao 
de R*/2,. 


Equivalência de sequências 


Vamos agora definir uma relação de equivalên- 
cia entre sequências reduzidas que permita, no 
conjunto cociente, incluir na mesma classe de 
equivalência as sequências que não difiram na 
natureza e ordenação dos estados, mas apenas 
na correspondência estabelecida com os elementos 
de R*/%, ou na constituição deste R+/9,. 

Poderíamos também referir directamente esta 
relação de equivalência às sequências 7, mas, 


RA RT 
havendo uma bijecção entre S, e Zp4, com fa- 
cilidade se passará de um a outro tipo de equi- 
valência. 


Definiremos então uma relação binária “* 
W Y . 
em =p; X 2p+ do seguinte modo : 


Sejam as aplicações 7; ,q;: r 4 tais que 
dq = (Ha ,Rº/90,ºSy) (20) 


= (Hi ,R'/255,5,) (21) 


q 
e 


Diz-se que 7, & 9; see só se existir um 
, Pt 

isomorfismo q de (R* /%s, <) sobre 
x 


(RÉ [ Pops É. ) tal que 
Cc =0 oq (22) 


Demonstra-se (Ver Anexo) que “* goza das 
propriedades reflexiva, simétrica e transitiva, pelo 
que é efectivamente uma relação de equivalência 
em 2r+, 

Paraleiamente diremos que duas sequências 
0 106 E - ai são equivalentes se se verificar 04 é 0%, 


é +! — . P--s R R 
em que é é uma relação binária em 5, a 5, + 
definida por 


R ” 
Dados 7,,7%€ 5, ”, diz-se que 5, & c; se 


e sóse 'y (74) 8* q (08) (23) 


e 


Dispensamo-nos de demonstrar que é é uma 
né nidos = Mc R + 
relação de equivalência em S, ',0o que resulta 
do facto de v ser bijectiva. 


Trajectórias 
Formemos o conjunto cociente 
T= Lor/S* (24) 


do conjunto das sequências reduzidas pela relação 
de equivalência acima introduzida. Cada elemento 
c:r é uma classe de sequências reduzidas, todas 
muútuamente equivalentes, isto é, redutíveis umas 
às outras através de isomorfismos entre as res- 
pectivas bases de referência reduzidas. 

Assim, em todas as sequências reduzidas 9* 
de uma mesma classe z (e também em todas as 
sequências 7 que lhes correspondem por ':) os 
estados de S ordenam-se pela mesma forma. 

Está completada a abstracção da base de refe- 
rência, desde que chamemos trajectória do sistema 
a cada classe de equivalência Z E 7. 
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Como classe de equivalência que é, cada trajec- 
tória z pode então ser representada por um 
qualquer dos seus elementos. 

Com esta passagem da sequência para a trajectó- 
ria, perdemos a possibilidade de estabelecer a 
simultaneidade entre um estado de S e outro 
qualquer fenómeno exterior a S e também perde- 
mos a noção do ritmo de evolução de 5. 

Conserva-se, porém, a noção de evolução 
cíclica que S possa ter registado. 

Note-se que as palavras «simultaneidade» 
e «ritmo» são usadas em sentido não necessária- 


mente temporal, mas apenas referido à base R* 
utilizada. 


5. COMPORTAMENTO DUM SISTEMA 
Regime livre 


Dissemos no nosso artigo já referido que, dados 
um sistema S e uma base de referência C, existe 
uma correspondência entre S e o conjunto das 
suas evoluções segundo C. 

Passando aos espaços imagem, diremos agora 
que existe, entre S,, imagem do sistema S em 

X Vi, e 2 + + Conjunto das sequências redu- 
je 
zidas, uma correspondência que, a cada estado 
observado E « S, associa uma ou mais sequências 
reduzidas 0*s : ZºR 4, tais que a imagem do pri- 
meiro elemento de R* /?ss é E — aquelas que S 
puder sofrer sem influências exógenas. 

A composição desta correspondência com a 
aplicação canónica de 3p + sobre 1,= 58 +/ 
/S*, gera uma correspondência entre 5: e q, que 
indica as trajectórias possíveis, em regime livre, 
com início em cada estado observado E. 


Regime forçado 


Da mesma forma, se o sistema S for sujeito à 
acção do meio exterior M, existirá uma corres- 
pondência entre S,X 342. + , e AS , que, a 
cada par (E, qu) : S,X 5Mp + associa uma ou 
mais sequências reduzidas 74 é: 2Sp +, iniciadas 
em E — aquelas que S puder sofrer quando M 
registar a sequência reduzida 04. 

Podemos da mesma forma, através da apli- 
cação canónica de XS, + sobre T«, passar a uma 
nova correspondência sobre o conjunto das tra- 
jectórias de S. Porém, isso tem aqui menos inte- 
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resse, na medida em que é necessário não per- 
der a noção de simultaneidade de estados de 
M edesS. 


NOTA FINAL 


Estamos agora em condições de poder passar 
à consideração dos sistemas como operadores no 
fluxo de informação, o que faremos em próximo 
artigo. 


ANEXO 


I) q 


q é uma relação de equivalência em R* 


a) A relação 9 é reflexiva 


Efectivamente, dado qualquer ieRt e qualquer 


GE SR” + temos que: 


(vi) [(iSk<)vli<ck<)— (=D) — 
—» (a (k) = c(1))] (25) 


o que, por (ro), prova que i VA lo. 
A relação 9 é simétrica 
Temos que demonstrar que, quaisquer que sejam 
ij:Rf e Gem”, se tem 
i É ini & i (26) 
Com efeito, temos então por hipótese, que 
(vk) [(i<k <DvG<Sk<Si) — o(k)=o(i))] (27) 
Ora é evidente que 
E SISN)vVISIS] (28) 
O termo da disjunção verificado depende de ser 
i<|jouj<i, podendo ser ambos verdadeiros se 
| x própria hipótese (27) impõe então que 
s (1) = a (i) (29) 
donde, ainda por (27) 


(vk) [i<k<)vi<Sk<SD)— (c(k)=o(])] (30) 


provando, por (10), que j Ji. 
c) A relação 3 é transitiva 


Temos que demonstrar que, quaisquer que sejam, 
,jkeR+ e c:SR+, se tem 


id,isjAk—igak (31) 


Com efeito, temos então, por hipótese 
(um [i<m<jpv(ji<m<i)— (a(m=e(i))] (32) 
(yum [(i<m<kyv(k<m< j)— (o(m) = c(j))] (33) 
Vimos em (29) que 
Pd — (c()=c(i)) (34) 
pelo que, de (32) e (33) vem 
(ym) [i<m<)vijga<am<)v(i<m<k)v 
v(k<m< j)-— (e (m) = (i))) (35) 
Ora, quaisquer que sejam i, |, k: R+, tem-se 
(i<m<k)—(i<m<jjvA<m<k) (36) 
ESm<)>(GSm<)vkCM<S]) (37) 
donde vem, por (35), que 
(ym)((i<m<k)vik<m<i)—(o(m)=o(i)jj)] (38) 
Isto prova que ij k, o que completa a demons- 
tração de que 9, é uma relação de equivalência em R*. 


2) Cada I:Rt/%, é um intervalo conexo em Rº- 


Sejal:Rt/J, e sejam dados ii: I, com i <i. 
Sendo I uma classe de equivalência segundo à e sendo 
i,i' «], temos que ig, À. 

Devido a (ro) pode então escrever-se 
(vk)[(i<k<il) — (o (k)=0(1))] (39) 

Consideremos agora | tal que 

ti E (40) 

Temos. por (39), que 

(vo [(i<k<p— (ick<Si)— (0(k) =ec(i))] 


(41) 
pelo que i 9,1. 


Mostrâmos assim que 
G<iDaA (te) — (vil) GD (42) 


o que prova que | é um intervalo conexo (fechado ou 
aberto) em RH. 


3) 88 é uma relação de equivalência em Zu+ 
a) A relação s* é reflexiva 


Efectivamente, dado qualquer 5º: > 4, basta tomar 
a aplicação identidade « de R?/2, sobre R*/9 a 
qual é um isomorfismo para a relação de ordem 
< , para provar que 


0” = e" o q (43) 


e portanto que 5º 8º 0º 


6 


b) A relação S* é simétrica 

Temos de demonstrar que, quaisquer que sejam 
* e 
dg 3 0B E 224 se tem 


* ek * * ex Nr 
Ts Õ 56 —— do! po Ta (44) 


Existirá então, pur hipótese, o isomorfismo « de- 


64 03 
(RF/%s, , <) sobre (RT/905 a <) tal que 


4 =9%609 (45) 


Mas, então, a correspondência inversa 2-1 de v é tam- 
“p 
bém um isomorfismo de (R*/%4,<) sobre 


Sa, 
(R$/%44, <<) e verifica, por az, 2e 9! serem 
aplicações 


oq (são?) oq”! 


= cão (7 oq) 


provando que 


c) Arelação S* é transitiva 


. % E 
Temos de demonstrar que, dados quaisquer q, 19%, 
7. : Zr+ , Se tem 


* S 


1% x 


o A 7; 6 1%. — To T+ (47) 


Ora existem, por hipótese, os isomorfismos y, de 


Ca Ga 
(RE / da, + <) sobre (RT / da, <)e % de (RE Jog, 
Ga - 


| = UA 
< ) sobre (R* /Y, ,<.) tais que 
o 


Ts mm 7 O e (48) 
* dr 


“ e e “ e 
Daqui resulta, por 7. , 2, e 7, serem aplicações, 
que 


pi = (55 O 2,) o 2% 


= q O) (2, O 2.) (50) 


Como ? 


a € ?, são isomorfismos, também 7, o 7, 


e] 


z 
é um isomorfismo de (R/9,, , <.) sobre 
6. 
(RT/9 8; + <. |, pelo que (50) prova que q, o 4. 
Fica assim demonstrado que “* é uma relação 
de equivalência em Les o 
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DA TEORIA DA FONTE MÓVEL NUM MEIO INFINITO 
A DUAS DIMENSÕES 


por  A.M.R. MARRECAS 
Bolseiro do IAC no NEEM 


RESUMO SYNOPSIS 

Generaliza-se ao regime transiente a solução da equa- The solution for two-dimensional heat flow due to a 
ção do calor no caso de uma fonte móvel com velocidade constant velocity moving heat source 18 extended to transient 
constante num meio infinito a duas dimensões. É usada a case, Laplace tranaform is used. 


transformada de Laplace, 


NOMENCLATURA 


T — Temperatura (origem 0) 

t — Tempo 

V — Velocidade da fonte 

k — Condutividade térmica 

« — Difusividade térmica 

OQ — Quantidade de calor libertada pela fonte na unidade de tempo 
(x,y) — Referencial móvel associado à fonte 

r2 iesna x + yº 

Jo ( ) — Função de Bessel de primeira espécie e ordem zero 

Yo ( ) — Função de Bessel de segunda espécie e ordem zero 

lo ( ) — Função de Bessel modificada de primeira espécie e ordem zero 
Ko ( ) — Função de Bessel modificada de segunda espécie e de ordem zero 
C,, Cs — Constantes de integração 


1— INTRODUÇÃO 


O processo de análise seguido em [1] é susceptível de, mediante algumas adaptações exigidas 
pelo novo modelo físico, conduzir ao resultado pretendido. A extensão ao regime transiente da 
solução correspondente ao regime estacionário que Rosenthal aponta em [2]. Esta a razão porque 
omitiremos observações e passos que nos pareçam redundantes com os inclusos em [1]. 


2 — MODELO TEÓRICO 
Difere do indicado em [1] em dois pontos: 


— (O meio onde se desenvolvem os fenómenos térmicos é um plano material isolado em ambas 
as faces; 


— Nos referenciais associados ao meio e à fonte omitem-se, agora os eixos dos ZZ. 


21 — Tratamento matemático 


A equação do calor não conterá derivadas em ordem aos ZZ, pelo que, escrita já no referencial 
móvel, reveste a forma 


OIT OT |V 9T 1 
a t (1) 


9T 
dx? 0 y? a dx 7 
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211 — Regime estacionário 


A integração da equação respectiva 


)? ê ) 
pts MME a (2) 
dx? dy * a dx 
conduz à solução [2] 
x 
24 
T= Q e ko (— ) 
27 k 2a (3) 
2.1.2 — Regime transiente 
A solução obter-se-á integrando 
9ºT Kimi Voo 3 wWtT (4) 
dx d y? q 0% a dt 
2.1.2.1 — Integração de (4) 
1) Condições fronteiras a respeitar 
Tso Para t=-0t, em todo o domínio (1.a) 
T-o Para t qualquer quando s — co (1.b) 
dT 
— 2 7 TKk Ee Q Para t>o qualquer quando r » o (1.c) 
dr 
2) Por proposição da primitiva 
V 
—— x 
T=e 2 fh(x,y,t) 
(4) Transforma-se em 
o? 9? 3 ) 
f 4 Eee (=) p= À df (5) 
dx? dy 2x a dt 
3) Aplicando a transformada de Laplace em ordem a t, (5) adquire a forma 
PE dE Erg S- 1 
—— — — (— | f=— E — + 
ds O dy 2a ; a T(o x,y) (6) 
em que 
- amem +00 — St 
E (x, y, Eros: | E(x, y, t) € dt 
o 
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TECNICA XIV 


E 4) Introduzindo a condição inicial (1. a) e praticando a mudança de coordenadas f (x, Y, 8) 
— f(r,s), (6) leva-nos a 


PE 1 Eid) + Sl f=o (7) 


5) Para integração de (7) mnotemos que se trata duma equação de Bessel modificada, admitindo 


Gt [yV(L) + +Sejrer [V(L) Ef | (8) 


6) Posto que lim lo VÁ =) + E | = 00 e tendo em vista (1.b), será C| = 0, dege- 


Tx, 


nerando (8) em 
é E Eos 2 A E ==: = 
Au MY. 5 
Ke [VD ++] (9) 


7) Uma análise das alíneas anteriores mostra que podemos, agora, escrever 


[run ce. 
V RA aj is | (10) 


8) Tomando em conta (1.c), modificada por aplicação da transformada de Laplace, ou seja 


E d T Q 
—2k7r dad quando r-— O 
e aplicando-a a (10) conclui-se ser 
Ca = Q 
2kxs 

9) V 

E DD E Tia Juve Ss | 

T = Ko j — — — 

2k= So |V 22 "ig ” (11) 


10) Na inversão da solução estará implicado, apenas, o conjunto 


ps 2 MIVYOo Ss 
» (r, 5) G Ko | V To “e = r (12) 
11) Determinação da função inversa 
die (ar sé à ” 
g (r, t) = = Jim, | g(r,s)e ds (13) 
Tv j-—+ 00 nicaa jB 


o contorno da Fig. 1, correspondente à expressão (13) será substituído pelo contorno da Fig. 2, resu- 
mindo-se as contribuições para g (r,t) a: | 
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Fig. 1 —(Y terá um valor suficientemente grande para dei- Fig. 2 — Assinaladas as singularidades da função inte- 
xar à esquerda todas as singularidades de g (r,s)Jest) granda 


— Resíduo da função integrada em relação ao polo S==0, o seu valor será 


2% 


— Integral estendido ao contorno BAEDC, a mudarça de variável 


2 it 
27 a 


permite transformá-lo na soma de dois integrais 


ao — q u?t aaa 
J [Jo (ru) —i Yo (ru)] e a a E a inô 
º 2 (aut + na) 
42 
"SO — q u?t E io 
À [Jo (ru) + i Yo (ru)] e a A e AM (16) 
2 (t+) 


Agrupando (14), (15) e (16), 
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V2 


V Ta né mai 
g (r,t)=ko pe ) ai f ue ço du (17) 
e” 0 PS . 
q u* + Za 
VZ Vv? DO 


2.1.2.2 — Conclusões 


1) (18) verifica a equação diferencial de partida (4) e satisfaz as condições fronteiras, dege- 
nerando em (3) quando t — 09, 

2) O integral contido em (15) deverá calcular-se por métodos numéricos. 

3) A solução tem interesse tecnológico, em particular no que respeita à soldadura de chapas 
finas por processos que envolvam uma fonte móvel. 
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ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL *(*) 


I—-Breve nota mensal 


As afluências ao conjunto do sistema revelaram-se, 
em Julho, sensivelmente inferiores aos valores médios, 
tendo totalizado apenas cerca de 127 GWh, « que cor- 
responde uma probabilidade de serem excedidas da 
ordem dos 85 "/, e um coeficiente de produtibilidade 
de 0,62. O carácter seco destas afluências ficou a de- 
ver-se, em parte, a um regime desfavorável de caudais 
transitados em Miranda, consequência de uma indispo- 
nibilidade total da central espanhola de Castro, por 
um período de 7 dias, durante o qual o caudal afluente 
a Mira foi praticamente nulo. 

As afluências acumuladas desde o início do ano hi- 
drológico têm uma probabilidade de serem excedidas 
da ordem dos 77 '/; e representam 0,82 do respectivo 
valor médio. 


NH — Elementos gerais (GWh) 


a) “Mensais Variação 
| 1967 | 1968 | 0, 

Produção hidráulica (Ph). 134,1) 360,1 1& 
Produção térmica (P+). 117] 120,x + + õB2 
Produção total (PT). ..... 156,8 | a O, 9| + 
Energia recebida de empresas 
não pertencentes ao ENC (Er) 1,5 UA e SO 
Exportações (Ex)... ...... 03 2,6 — 
a repsiç TEA é ss mis » é 1,2 9 + 58 

Saldo importador (S1). + 09— 05 esa 
C one. em bomb hidroeloé JC 5) (1) 0,0 0,0 0 
Produção para con- | . (1) 

sumos perman. (Pcp) (7) | 3851 444 +11, 
Produção para con- . | 

sumos não perman. (Penp) () , 431) 532,— N 
Total. o 45» E 40,61 + 5 
Coeficiente de hidraulicidade 0,73! 0,62 L 


NOTA: 

(1) O aumento percentual da produção para consumos permanen- 
tes, tendo em conta a incidência dos dominges e dias especiais é 
respectivamente de 0,5 e 11,1 So» 


[HI — Diagramas de cerga dos dias cerrcterísticos 


em 


ma 4,* feira: 


NM 1968 

Produção hidráulica (Ph) MWh | 16187 [13066 
Produção térmica (Pr) MWh 1) 8 230) 
Produção tota! (Pr) MWh 1615716296 
Trocas com | Export. (Ex) MWh 0 100 
Espanha | Import. (1) MWh/ O 20 
“Cousum em bomb. hidroel. (Cb) uun O Ú 


(ep) MWh [18651 [14695 


Prod, para cons. perm, 
2 441 1521 


Prod. para cons. não perm, (Eenp) MWh 


TOTAL Pi +(IEx) M MWh | 16142 [16216 - 
Potênciamás. MW. 85] eso 
Potência min. MW 488 436 
I-E 
EF PTH OD | Ueiliz, da ponta horas 18,3 18,3 
z É | Factor de carga 0,76 0,76 
Es [Potência max. MW. T7S | std 
Ea p Potência min. MW. 331 304 
io “p Utiliz, da ponta horas 17,0 18,0 
Factor de earea 0,71 n,75 


[V — Energia «rmazera1a mas principais aliufeires 


No fim do mes 
Albafeiras: 


GWh 0,0 (1) 
Alto Rabagão . . 2... «| 4526 46,5 
Paradela . 155,2 64, 
Venda Nova +. +... «0 120,6 94,2 
Salamonde . . .. «o. sl 26,4 45,1 
Cariçe: RES NE , 28,4 55,8 
VE 0a. Eai 51,1 53,5 
Cabril .. Ema 251,8 14,2 
Castelo do Bode... 105,7 (4,8 
Guilhofrei . SE E | 6,2 Ta 
Lagoa Comprida 39,4 (2) 86,9 
Santa Luzia +... cv. 06,9 91,7 
FrICSNE: É ug Gus s%ê Z Na | 6,0 
Povoa E xe é 11,2(3) 89,4 
Total | COM A, Rabag são o 1 313,4 61,9 
sem À, Rabagão 860,1 4,9 
NOTAS 


(1) Coeficiente de enchimento em energia das albufeiras, 


(2) Inclui 3,2 GWh armazenados em Vale do Rossim no início do 
mês e 28 GWh no fim do mês. 


(4) Inclui 1,0 GWh armazenados no açude do Poio no início do mês 
e 1,7 GWh no fim do mês. 


(+) Elementos extratos das estatusticas mensais do Repar- 
tdor Nuciunal de Curgas (RN. C,). As produções e 
os consumos das empresas do R. N. C. representam 
cerca de 940/, dos totais do Pais, 
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dist Belarte 


AEG-LUSITANA DE ELECTRICIDADE 


Qualquer destes nove-tipos 
de disjuntores é mais rápido do que 
o curto-circuito 


Os disjuntores AEG-RAPID são disjun- 
tores limitadores de corrente, A sua 
rapidez de actuação e tal, que qualquer 
curto-circuito é eliminado logo à nas- 
cença 


Os disjuntores AEG-RAPID limitam a 
grandeza e a duração duma corrente de 
curto-circuito à valores inofensivos Com 
cstes disjuntores, a instalação iica so- 
mente submetida a 20%, de solicitação 
dinâmica, ao contrário do que sucedia 
com os disjuntores de construção mais 
antiga. Assim, deixam de venrficar-se as 
usuais quebras de apoios e deformações 
dos barramentos. Há somente que con: 
tar com 10% de sobrecarga termica. 
eliminando-se assim incêndios em cabos 
e celas. Por outro lado, não surgem dete- 
rorações resultantes de arcos voltaicos 
migratórios, na medida em que o curto- 
circuito fica limitado ao local onde se 
onginou. 


Como vantagem consequente teremos 
Produção processual sem interferências 
e sem avarias secundárias. O disjuntor 
AEG-RAPID é fabricado para a gama de 
I6A a 20004, 500V 3” e um poder de 
corte até IO00kA 


Aproveite V a oportunidade, elabore os 
projectos do futuro, recorrendo a disjun- 
tores de avançada concepção tecnica 


A AEG estã sempre, com prazer, ao v/ 
dispor para quaisquer esclarecimentos 


em toda a parte 
aparelhagem de corte 


AEG* 


um sinal de progresso. segurança e confiança 
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ESTAMPARIAS TINTURARIAS 
LAVANDARIAS 


PARA OS VOSSOS TECIDOS TÉCNICOS: 
TRANSPORTADORES, MANCHONS, LAPPING, WOLFRIES, 
FLANELAS, CANEVAS, ETC,, 


EM FIBRAS NATURAIS OU SINTÉTICAS 
CONSULTEM: 


FANAF EL 
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O EFEITO DAS VARIAÇÕES TÉRMICAS SINUSOIDAIS 
NAS ESTRUTURAS LAMINARES DE COMPORTAMENTO 
ELÁSTICO LINEAR 


SUMÁRIO 


Us estados de tensão criados numa peça laminar de betão 
pelas oscilações térmicas ambientes podem ser analisados 
por métodos numéricos, se se considerarem aceitaveis deter- 
minadas hipóteses relativas à transmissão do calor e às 
propriedades termo-elásticas do betão. Depois de se formu- 
larem as bases para esta análise, desenvolve-se um metodo de 
cáleulo numérico em termos das componentes dos desloca- 
mentos de pontos discretos duma casca de geometria arbi- 
trária. 

O metodo é exemplificado na análise do estado de tensão 
duma barragem casca submetida às oscilações anuais de 
temperatura, 


NOTAÇÕES 


X, y,Z, X, Y, Z — coordenadas cartesianas orto- 
gonais 

é, ',E, - coordenadas curvilíneas ortognais 
l,m,n  cossenos directores duma normal à su- 
perfície média 

uv, w, U, V, W — componentes do desloca- 
mento 

t — variável tempo 

9 — temperatura ('C) 

E — extensão 

7 — distorção 

Z — tensão tangencial 

E — módulo de elasticidade longitudinal 

G — módulo de elasticidade transversal. 

v — coeficiente de Poisson 

a — coeficiente de dilatação térmica 

N — esforço normal 

T — esforço tangencial 
M — momento flector 
4 — momento de torção 
Q — esforço transverso 


por unidade de compri- 
mento da peça linear 


* Engenheiro Civil, LNEC, 


por JJ. F. LOBO FIALHO * 
M. MACHADO GONÇALVES ** 


SYNOPSIS 


The state of stress in a concrete shell due to the surroun- 
ding temperature changes can be analyzed by numerical 
methods if some hypotheses concerning heat transmissions, 
and thermo-elastic characteristics of the concrete are accepted 

In the present work the basis for this numerical analy- 
sis is developed for shells with arbitrary shapes. The propo- 
sed method is exemplified in study of stress in a shell dam 
due to the annual temperature variations. 


1 — GENERALIDADES 


O efeito das variações térmicas nas estruturas 
laminares de betão ou betão armado utilizadas 
na construção (coberturas em casca, lajes, abóba- 
das delgadas, etc.) pode ser feito em primeira 
aproximação, admitindo que são sinusoidais as 
oscilações da temperatura ambiente. Conforme a 
espessura da estrutura laminar, assim a predo- 
minância da onda sinusoidal diurna ou da onda 
anual no estado de tensão gerado pelos gradi- 
entes térmicos. 

A fixação da onda térmica predominante, isto é, 
do seu período e da sua amplitude, depende de 
numerosos parâmetros geométricos e de exposi- 
ção da obra, bem como de condições meteoroló- 
gicas ambientes, que não serão analisados no 
presente trabalho. Este visa fundamentalmente o 
estudo do estado de tensão gerado numa peça 
laminar de geometria arbitrária, supostas conhe- 


cidas as oscilações sinusoidais da tempera- 
tura ambiente na proximidade das suas super- 


fícies. 


** Engenheiro da Universidade Católica do Rio de Janeiro, Bolseiro da Fundação C. Gulbenkian. 
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Para realizar este estudo consideraram-se as 
seguintes hipóteses: 


a) A peça laminar é homogénea e simples- 
mente conexa. 

b) O material de que se constitui é isotrópico 
e linearmente elástico, mantendo estas 
propriedades, bem como as constantes elás- 
ticas dentro da amplitude de oscilação consi- 
derada. 

c) O problema termoelástico é do tipo “quase- 
-estático”, isto é, o período da oscilação 
térmica é tão longo que se podem despresar 
os fenómenos de inércia e os termos de 
acoplamento termo-mecânico [1]. 


As hipóteses (a) e (b) são admissíveis, visto 
que o betão e o betão armado, dentro das ampli- 
tudes de oscilação térmica ambiente se podem 
considerar como um comportamento hookeano. 

A hipótese (c) tem sido aceite por numerosos 
autores quando as variações térmicas são lentas 
comparadas com a velocidade de transmissão do 
calor através da peça laminar. Ela permite simpli- 
ficações importantes nas equações de equilíbrio 
tornando possível um tratamento numérico do 
problema. 


2 —- ESTADOS DE TENSÃO E DEFORMAÇÃO 
TERMOELÁSTICOS 


Considere-se uma peça estrutural a três dimen- 
sões referida a um sistema de eixos ortogonais 
cartesianos OXYZ suposta no estado neutro à 
temperatura constante “o. Por estado neutro 
entende-se aquele em que os elementos de volu- 
me estão todos indeformados e sem tensões 
iniciais [2]. Suponhamos que esta peça é subme- 
tida num dado instante a uma variação de tempe- 
raturas 0 ==) —0, uniforme em todos os seus 
pontos. Se a peça não tiver ligações que impeçam 
a livre dilatação de qualquer dos seus elementos 
de volume, isto é, se não houver fenómenos de 
coacção e além disso a variação de temperatura 
se fizer em regime “quase-estático”, ter-se-á um 
estado de deformação total 


7 ——— 0, a 0 (1) 
a que correspondem as expressões escalares 
(1 a) 


14 


a do a O! b) 
As tensões normais e tangenciais em qualquer 
ponto serão portanto nulas visto que todos os 
elementos de volume se podem dilatar livre- 
mente sem efeitos de inércia. Se, pelo contrário, 
como será em geral o caso, for criado um estado 
de coacção às dilatações térmicas devido a uma 
distribuição não linear de temperaturas ao longo 
de uma secção da peça ou a uma ligação de coacção, 
aparecerá na peça um estado de tensão de coacção 


nx não nulo e um estado de deformação termo-elástico 


E” diferente do dado por (1). O estado de tensão 
de coacção é resultante dum estado de defor- 


o Q 
aos 1a 
o =] 
<a É a 
| ——ee e e em L-—-— 4 
K € 
] E [3 
MM 
> aê Si 
O o] 
sa q — 
Usus 
n5 Au 
To E 
E-—oe---e--s 
| 
O 
Cc a 
E meo 
Le 
Fig. 1 


mação de coacção CS sida (Fig. 1) a que 


ij 
corresponde a expressão tensorial 


EC = 68.0 n 
ii tj Ria (2) 


Esta expressãs pode desdobrar-se nas 6 relações 
escalares 


Es E (2 a) 
Cota 
Cima io 
Pe Lug 
Eee: dp: É 
Su 
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Se o material da peça estrutural se mantiver 
isotrópicamente elástico conservando as respec- 
tivas constantes, ter-se-á: 


É ie c 
“ii ami ui (3) 


a que correspondem as expressões: 


E 
ge 


(1-2) (14) (1 did | 


E 


= a) qo 094 ++] Ga 


“O (1.2) (149) [à rãs, | 


Csy = G Vas 
E =6 1h (3 b) 
sz = 6 pl 


Substituindo os valores de :Z pelas suas expres- 
sões (2 a) e (2 b) ter-se-á: 


co. E toe O É e 
+ (1-2 v) il ra TR y E) 


= (1 +») a0 |=a? “Eat 
1-2 » 


E f 
al =— — 1-2) c 8 v (e 0) — 
d pra! niá added 
(4a) 
-049) 40 |=oj-Eos 
1-2 v 
E = [Genebra 9 pb 
* (1-2) (1 +) á » 
-0+907 |=op Ee 
1-2 » 
Li] 
9 ay 
cm a b) 
Ze femea 4 
ye Cy 


Sendo cj o estado da tensão que produz na 
peça o estado de deformação termo-elástico. 
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O estado de tensão de coacção deverá satis- 
fazer as aquações de equilíbrio indefinido quando 
nulas as forças de massa, isto é: 


d 0x d 4 ) 4 

A od 

d 0y d Cas d Ta a 

» o PS (5) 
Ô oz Ô Txz d Tz o 

dz + dx SP dy ? 


e às condições de equilíbrio na fronteira quando 
nulas as forças de superfície: 


log mz +tncç=o 


Zso + m 9€ 


; ENTE=O (6) 


tm +tnos =0 


Se substituirmos em (5) e (6) os valores das 
componentes de c;j pelas suas expressões (4 a) 
e (4 b) obteremos as equações de equilíbrio inde- 
finido : 


0% OM PM E OT, 
ox dy dz (1-2) dx 
o 09 dt Rr Ai MM 
Epotrgo a Pg uq ff 
Ty dx dz (1-2v) dy 
f — 
05 o, A a - 498 ca 
dz dx dy (1-27) dz 
e as equações de equilíbrio na fronteira 
E Pe 
A Etta nto 0=0 
lz) +m? + A x0=0 (8) 
Ny “m y ni Zy (1— 27) RR ic 
Iz)+mz! + ed men ii a 0 =0 
xz yz z (1-2) 


A análise das expressões (7) e (8) permite 
enunciar a seguinte analogia atribuída a DuHa- 
MEL [1]. 

O estado de deformação termo-elástico é idêntico ao 
que se observaria se a peça fosse solicitada por forças de 
massa 

E d 6 
x 
o Xi 


X, = 


i fi (Tm 27): (9) 


e forças de superfície 


E 


——— (10) 
1—2v) 


lag 


O estado de tensão de coacção obtém-se das 
expressões (4a) e (4b) que têm na notação 
tensorial a forma: 


E 47 (11) 


3 — CASO PARTICULAR DAS ESTRUTURAS 
LAMINARES 


A analogia anteriormente exposta pode ser 
generalizada ao estudo das variações térmicas 
das estruturas laminares. As respectivas expres- 
sões devem porém ser formuladas nas variáveis 
habituais da teoria das cascas. 

As equações de equilíbrio indefinido em ter- 
mos dos 8 esforços de casca, N;, N,, Ti, M:» 
M,, VW, ,/Q:e Q,, segundo as coordenadas É, 7 
e É têm a seguinte forma [8]: 


DN; |, DT, , N—N,; 2T, 4 
DE Dn R; R* 
Q Q, 
S ata: + f; = O 
Ry R 7» 
DN, DT, i N—N, 2T, E 
Dn E Rº R; 
Q, Q; 
- — em —+f,=0 (12) 
RE? 
om é PÕE 
OQ Q, 
Ei o Ei z ci É; dio 
E 


Sendo 1R;, IR” e 1/R$ componentes do ten- 


sor de curvatura nas direcções é e 1; Rº eR, 
os raios de curvatura geodésica no ponto consi- 
derado e [f:, f,, f | o vector solicitação refe- 
rido à unidade de área da superfície média; 
D/D: e D/Dn são derivações totais (substan- 
ciais) ao longo das abcissas curvilíneas ; e n 


traçadas sobre uma rede de curvas paramétricas 
ortogonais existentes na superfície média da 
casca. 

Entre os esforços de flexão M; e M, , de tor- 
ção ol, e os esforços transversos (OQ: e Q, exis- 
tem as seguintes relações de equilíbrio : 


DM: | DM,  2M, M,—M: ES ei 
D& D n KR: R$ e 
(13) 
Dn DE Rê Rº , 


Substituindo os valores de O; e Q, dados 
pelas expressões (13) nas equações (12) obtêm-se 
3 equações diferenciais relacionando apenas os 
esforços NE, N,, ME, M,, Tr, € Mk. 

Entre estes esforços generalizados e as corres- 
pondentes deformações generalizadas* da super- 
fície média, & ,€,, kz , ky, 7x, TE e T,, existem 
as relações 


NE =K, (gw +v6,) 
N,=K,(º, +v%) 
Ten = Kg yr 

M,=K,(k: + vk,) 
M,=K,(k,+vk:) 


(14) 


OW, meto K, (42 aan K, Cy 


sendo K,, K,, Ke K,, os módulos de rigidez 
axial, tangencial, de flexão e de torção, dados 
respectivamente pelas expressões: 


q Es (15) 


(*) cE + E, — extensões da superfície média segundo 
as direcções £ e 7. 
Y:n — distorção da superfície média nas direc- 
ções 5 e 1. 
kz ,ky — variações de curvatura da superfície mé- 
dia segundo & er. 
T: , Ty — variações de torção segundo £ er. 


- 4 
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3 
K, = te 
12 (1—v2) 

3 

EE E e 
24 (1 ++) 


Postas estas questões de base da formulação 
geral da teoria das cascas, analisemos o que se 
passa na estrutura laminar quando esta sofre 
uma variação de temperatura 9 (£,7) num ponto 
genérico da superfície média com uma variação 
(linear) 9 (£, n) entre os dois pontos situados 
nos paramentos (Fig. 2). 


e. —— 


Fig. 2 


Em virtude desta variação de temperatura, as 
deformações de coacção em cada ponto da super- 
fície média da casca são: 


6; = 8, E) 


4,=8 —ab (16a) 
Í| 
Ent 
c / E] 
kp=kp— &- 
e 
Et. 4 
ke =kç— a — (16b) 
e 
C f 


onde Sº , y', kºe zº são as componentes do 
estado de deformação termo-elástico. 

Note-se que as expressões (16 a) e (16 b) 
correspondem às expressões (2 a) e (2 b) da for- 
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mulação elástica tri-dimensional em coordenadas 
cartesianas ortogonais. 

Substituindo (16a) e (16 b) nas expressões 
(14), obtém-se: 


NE=N;—K, (1+7)a 0 


Ni=N!—K(1+2)20 (7a) 
TE T; 
! x 6 
M; = E = K, (+) ã 
re ab , 
MM. — K (+?) (17b) 
e 


sendo Nº, Te, Mº,e Nº os esforços de coacção 


eN',T',M e ou os esforços que produzem 
o estado de deformações termo-elástico. Note-se 
que as expressões (17) correspondem às expres- 
sões (11) da formulação cartesiana. Substituindo 
as expressões (17) nas equações de equilíbrio (12) 
depois de eliminados os valores de O; e Q, a 
partir de (13) e fazendo-se ainda f; = f, =f, =, 
obtêm-se três equações de equilíbrio em tudo 
análogas às equações de equilíbrio que se obteriam 
se as componentes da solicitação actuante fossem: 


e o K,º 
ionrge) PRO, 1 Pt) 
S D & R; D É 
K, 
la 
Rºé Dm 


(18 a) 


17 


| 0 0 
tfi=(+vya Kat ka Ji 
RR. 
mas E q 
E K, 9 
1 Pl) O 
O e Z -t- EEE USE o 
RP EM 
0 K, 
D(——) D(——) 
- RSS, 
D E D nº 


Estas forças fictícias podem considerar-se como 
correspondendo nas cascas às forças de massa (9) 
da analogia de DumHaMEL nos equilíbrios elás- 
ticos tri-dimensionais. 

No caso da estrutura laminar possuir um 
bordo livre normal à direcção arbitrária no, de- 
verão observar-se nele as seguintes condições de 
fronteira : 


as quais, substituídas em (17a) e (17b) corres- 
pondem à aplicação de forças fictícias no bordo 
livre: 


o = (1+7) ak (18c) 
e a momentos fictícios no bordo livre: 
hi K y 
Mo = (1+%)a— (18d) 


e 


que podem considerar-se corresponder nas cascas 
as forças de superfície (10) da analogia de 
DUHAMEL, 

Note-se que em geral se pode considerar no 


bordo livre 1 ==0, donde Mio = 0. 


4 — DETERMINAÇÃO DAS TEMPERATURAS 
Ta 


A formulação desenvolvida no parágrafo an- 
terior permite calcular em cada ponto as com- 


y 5 6 à gi rá Esá 
ponentes f. , o ef, da solicitação fictícia que 


produz, juntamente com as forças de bordo, na 
peça laminar, deformações idênticas às defor- 
mações termo-elásticas. As expressões (18) mos- 
tram que este cálculo só pode ser feito desde 
que se conheça a distribuição das temperaturas 


0 e 6 em cada ponto da peça laminar. 

Sendo a variação de temperatura ambiente 
sinusoidal de pulsação w prova-se [3] que tam- 
bém o serão as temperaturas 5 e 1, possuindo 
estas a mesma pulsação. 

Suponha-se que a peça laminar é uma placa 
indefinida de espessura constante e . Sejam: 


Om == 0 m + Am cosnt 
(19) 


=0;+ A; cost 


as variações de temperatura ambiente do lado 
da face m e da face j (fig. 3), as quais têm 
respectivamente as amplitudes Am e A;, oscilando 
em torno das temperaturas médias Om e 2. 


Integrando a equação de Fourier da trans- 
missão de calor aplicada à placa indefinida, 
supostas as suas faces submetidas às variações 
térmicas ambientes, expressas por (19), RITTER 


obteve [4| para um ponto interior de abcissa x: 


HJ 


0=NV9? + 4? cos | ot —are tg (5) dm 


no. no 
+ ( 57 Bm ) x (20) 
e 
sendo : 
2 == Am (cosux shux + B' cosvx shux + 
+ C' senyx chux) (21) 
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y == Am (-senux shux — B' senux chux + 
“+ C' cosux shux) 


n cosue — chue 


em que 


'—shpe che — costpe 


n shue — cosue 


. = 22 b 
ci che — cos?ve ae) 
À; 
= e 23) 
n An ( 


e » uma constante relacionada com a massa 
específica y, o calor específico C, a condutibili- 
dade térmica ) e a pulsação w, pela expressão: 


(24) 


Vê-se assim que a temperatura num ponto do 
interior da placa varia sinusoidalmente com um 
período T=2 n/nw pj torno duma temperatura 

rá 
média 0 = Om + je 


— E x com uma ampli- 


tude A==V 9? + 42, apresentando, em relação 
às oscilações da temperatura ambiente supostas 
em fase de um e outro lado da placa, um desfa- 


“ 
samento igual a arct. tg = ). 


O facto da distribuição das temperaturas ao 
longo duma normal à placa se fazer segundo 
uma lei não linear (fig.) 4) que varia de instante 


8,-temp? de referência 


Fig. 4 
para instante levou RITTER ao estabelecimento de 


dois artifícios para o calculo das temperaturas : 
equivalentes ge 9: 
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a) Procurou uma distribuição linear de tempe- 
raturas ao longo da espessura da placa 
de efeitos equivalentes aos da distribui- 
ção real, isto é, produzindo a mesma trans- 
lação e a mesma rotação duma normal 
à placa. 

b) Considerou, dentro dos diagramas lineares, 
apenas o mais desfavorável, isto é, o que 
correspondia à maior translação e rota- 
ção. 


O artifício (a) permite a determinação da tem- 
peratura média 6 do diagrama linear equivalente : 


1 [“'e/2 
) =— I) ndx (25) 
e 
— e/2 
e a sua variação 9: 
= 12 e/2 
= + — ) 9 xdx (26) 
é 
—e/2 


Substituindo em (25) e (26) os valores de 9 
dados por (26) e calculados os valores extremos 


de 9 e 9 a partir da anulação das respectivas 
derivadas relativamente à variável tempo, tem-se: 


o E 
air: = — VC f xá po x dgtA 


—e/2 — e/2 
9º + 9€ 
+——S (27) 
À 
o “e/2 o 
* extr = +84] Poxde [Cond] bx dx)? + 
raca ama 
0º — 6 
+ dE 28 
= (28) 


sendo, de acordo com (21): 


"e/2 Am(1+n) shze-+senpe 

J odk=T"—— >—————— (29) 
2p chue+ cosue 

— €/2 


e/2 Pa (1 + n) 
| Es E e cm 
ú 2u chge +'cospe 


shue—senpe 
(30) 


— e/2 
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e An(I—nle sh “ue—senue 


| md muco CU 
— €/2 
[Oy rm Am(1— n)e | sh e+-sen ve E 
44 chze—-cosve pe 
— e/2 
(32) 
Note-se que embora os valores de or. 


e Oextr. não ocorram em geral simultâneamente, 
esta hipótese poderá ser sempre considerada por 
se situar do lado da segurança. O artifício desen- 
volvido para as placas pode considerar-se ex- 
tensivo às peças laminares de pequena curva- 
tura, isto é, às cascas de espessura pequena quando 
comparada com a sua curvatura média. É o que 
vai fazer-se no parágrafo seguinte para o cálculo 
do vector das forças fictícias. 


5 — CÁLCULO NUMÉRICO 


Fixadaa lei de variação da temperatura am- 
biente nas duas faces da casca, a formulação 
estabelecida no parágrafo anterior permite a 


determinação das temperaturas 0 e 9 mais des- 
favoráveis que deverão servir de base ao cálculo 


das forças fictícias É: ; E , £ e fo - Às ex- 


pressões (18) podem ser traduzidas em diferen- 
ças finitas, o que significa que a determinação 


das temperaturas 0 e O se fará apenas para uma 
rede discreta de pontos. Sejam 9 (0,0) e 9 (0,0) 


respectivamente os valores de 0 e 0 num ponto 
central «== 0, $ = 0 tomado como referência 
numa rede ortogonal ax É traçada sobre a 
superfície média da casca. Representemos os 
pontos situados à frente e atrás de (0,0) sobre 
«== cte., respectivamente por (2,0) e (— 2,0) e 
sobre É == cte. por (0,2) e (0, — 2), sendo as dis- 
tâncias respectivamente (E> mm Ea) e (15 — no2) : 
Adoptando análoga notação para os pontos vizi- 
nhos as expressões (18a) e (18c) tomam em dife- 
renças finitas a seguinte forma: 


(Ka 9 ) 2,0 — (Ka 0). 2,0 


E (0,0)= (1-7) à | | 
i Li E—6,) 
Kr5  Ky0 
4 1 l e Fire 
= (E — E) 
t 
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RE) (Re) 


e e o-s 


A 
= 


(15 — nº 9) 


(Ka DoglKa D 94 


E. (0,9) =(1+92) 


(15 e 12.4) 
1 Kç8 K. 9 
Ru 
v 0,2 e “0-2 
= o o im aeniás 
Rº (12 — n-9) 
Ke 9 o K,0 à = 
a U e | 20 y e | 20 
Rº VR A 
, (Ka 0100 (Ka 9)0.0 
vi = ( DD +t-———— — 
fr (0,0) = 11 + v) « | R? Rº 


[ K$ 0) 2 Kç0 , 
EPE! 


e E A me (4 
(& =) —t9 


(6 (8 


e 0,0 


E ts E = E 


Ke 9 ú K 0 
eum U jp 


(ng — o) (rg — 123) 


rea 
e bo e 0, — 4 


+ = Sa 


EP) 


e 


+. 4 ço r0 = + 
R' (5, e 5.9) 
( Ki ) nã, te dy ne” 
dy 1 e mu. e 0,—* 
Rº (n2 — no) 


fio = (1 + v) a K,º 


TECNICA Nº 381 


CONSTRUCOES 
TECNICAS, LDA. 


PRAÇA DO MUNICÍPIO, 13-3º TELEF. 366506 LISBOA 2 
DELEGAÇÕES: CABO VERDE - ANGOLA - MOÇAMBIQUE 


EMPRESA TERMOELÉCTRICA PORTUGUESA 


Ensecadeiras para execução da tomada de água da Central do Carregado. 


TÉCNICA — X VII 


CONSTRUCOES 
TECNICAS, LDA. 


PRAÇA DO MUNICÍPIO, 13-3* TELEF. 366506 LISBOA 2 
DELEGAÇÕES: CABO VERDE - ANGOLA - MOÇAMBIQUE 


mm 
UU 


4 Ta 
e 
-—- 
ms 
Sa tt dr 
gr: tas 
a = E 
em Rê 
== ES 3232" Eos - 
e - == === 25 28 == 58 | | Bda 
=== 
o mu 
ES ES FS ES 
== 33 35 33 E 
e T— = = 
À | 
= a “ 
= i | 
E ! 
: rm s E 
. ” E 


— 
GEES | 
a Eu “TR 


dm. a À 
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E PORTUARIAS E VIADUTOS 


CONSTRUCOES CONSTRUCOES 
CIVIS INDUSTRIAIS 


Como se disse no $S 3, as expressões das equa- 
ções de equilíbrio (12) e (13) podem ser reunidas 
em 3 únicas equações diferenciais que relacio- 
nam os esforços de casca No. N,M,, M,, 
T, = Tu: € Ha == e, . Por outro lado as ex- 
pressões (14) permitem substituir naquelas equa- 
ções de equilíbrio os esforços pelas deformações 
generalizadas es o k. , k, rt Tel, OU 
ainda estas pelas três componentes u, v e w do 
deslocamento, através das expressões [8]: 


Du v w 

a = = 
D & R, R » 

pi ais (34) 
Dn R: Ç 


A we D(2)+ 
DaD; R; D é D&t 


D'w 1 Dw D u 
e tor Do Fo (=) 
D: Dr R: Dr; + Dn ( + 


v D v ) u 
ps: e Ss) uBnfb 
CRE O Di R? RR; 


As três equações diferenciais, de equilíbrio 
assim obtidas em termos das três componentes 
do deslocamento permitem a resolução analítica 
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do problema, desde que se considerem como 
incógnitas os deslocamentos. É fácil verificar que 
o sistema, das três equações diferenciais de equi- 
líbrio é de 8.2 ordem, sendo duas das equações 
de 2.2 ordem e a 3.2 correspondente ao equili- 
brio segundo a normal à casca, de 4.2 ordem. 

Aquelas equações podem também ser expressas 
em diferenças finitas, obtendo-se um sistema de 
3 N equações a 3 N incógnitas, em que Né o 
número de nós discretizados sobre a superfície 
média da casca. 

Utilizando a notação matricial pode dar-se ao 
sistema de equações obtido a seguinte forma: 


KU + fº =0 (35) 
sendo K uma matriz quadrada de dimensões 
3Nx>3N,se o vector das três componentes do 


deslocamento de nós considerados u = ! (u,, Vj, 
Wi) (UG, Vas Wa), -.. (Og Vo Wa)] e £º o vec- 


tor das três componentes Era f x e E: das força, 


fictícias correspondentes aos mesmos nós, incluí- 
das as de bordo livre. Vê-se assim que desde 
que se disponha da matriz K, formada por 
exemplo para o estudo de outra solicitação funda- 
mental actuando na peça laminar, o estudo do 
efeito das variações térmicas é facilmente reali- 
zado. Com efeito, basta calcular as componentes 


do vector da solicitação fictícia fº correspondente 
às variações térmicas e efectuar a operação 


u=—K”!gº (36) 


ou seja, multiplicar o vector — fº pela matriz 
inversa de K. Conhecidos os deslocamentos nos 
nós, obtém-se pelas expressões (34) e (14) os 
valores dos esforços. Note-se porém que os 
esforços calculados são os esforços termo-elásticos 
actuantes na peça laminar, correspondendo por- 
tanto aos assinalados com o expoente / nas expres- 
sões (17). Para se obterem os esforços de coacção, 
isto é, os que produzem o estado de tensão na 
casca, há que aplicar aquelas expressões a cada 
um dos nós considerados. Note-se que nesta fase 
há que aplicar ao bordo livre esforços que anu- 


lem os esforços N/ impostos ao contorno. Conhe- 


cidos os valores dos esforços de coacção, obtêm-se 
imediatamente as tensões instaladas na su- 


perfície da peça laminar devido ao efeito tér- 
mico : 


Nº 6 M: 
; É ca À 
N; 6M; 
q = mes (37) 
e — e 
T: 6 //8 
Ens = sp - 
e — €2 


6 — EXEMPLO 


O método de cálculo exposto vai ser aplicado ao 
estudo do efeito das variações térmicas ambientes 
sobre uma barragem casca. As características 
geométricas da barragem escolhida são as defi- 
nidas no trabalho anterior dum dos autores [8]. 


Fig. 5 


A discretização realizada sobre a superfície média 
da barragem casca definiu 75 nós pertencentes 
a uma rede ortogonal de curvas, paramétricas, 
englobando como curvas da mesma família o 
coroamento e a linha da fundação e como curvas 
da família que lhe é ortogonal a que existe no 
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plano vertical de simetria da barragem e a dos 
encontros. 

A matriz rigidez já formada para o estudo do 
efeito da pressão hidrostática [8] foi de novo 
utilizada para o estudo do efeito das variações 
térmicas. Na figura 5 apresenta-se o valor das 
temperaturas médias 9º e as amplitudes A da 
oscilação anual ao longo dos paramentos de 
montante e de jusante. Supõe-se que a tempera- 


tura média 9º é constante em ambos os paramentos 
+ 15ºC, havendo amplitudes máximas de 
oscilação anual de + 15 ºC no coroamento. Estas 


amplitudes mantêm-se constantes do lado de 
jusante, mas decrescem linearmente com a pro- 
fundidade, a montante, até anular-se no fundo da 
albufeira. 

Estes diagramas esquematizam os resultados 
de observações efectuadas em numerosas barra- 
gens construídas em climas temperados [5]. A 
partir deles obtiveram-se, pela formulação esta- 
belecida anteriormente, para os diversos nós da 
malha discreta, as forças fictícias indicadas no 
quadro I. Neste quadro apresentam-se também 


os valores das temperaturas 1 e 1, bem como 
ourtos valores que serviram de base ao cál- 
culo daquelas forças. É importante notar que os 


valores de  figurados no quadro estão referidos 
à temperatura de injecção das juntas de beto- 


nagem, que foi tomada como sendo /, = + 10 SE. 


A partir do vector de forças fictícias realizou-se 
no computador ELLIOTT 803 do L.N.E.€. 
a operação matricial expressa por (36) obten- 
do-se as componentes do vector deslocamento 
apresentadas nas figuras 6, 7 e 8. 

Utilizando programas já elaborados, calcularam- 
-se os esforços de casca correspondentes aos 
deslocamentos termo-elásticos. Estes esforços 
fictícios têm que ser corrigidos pelas expressões 
(17a) e (17b) para se obterem os que originam 
o estado de tensão de coacção. Também nesta 
fase foi desfeito o artifício utilizado para traduzir 
o comportamento do bordo livre, aplicando for- 
ças iguais e oposta àquelas que se impuseram 
aquela fronteira. Na figura 9 apresentam-se os 
valores dos esforços N.,N ,M. e M, ao longo 


do perfil central (2 == 0). Na mesma figura indi- 
ca-se a tracejado o andamento dos esforços termo- 


à q 1) f 5 ú. , a E 
elásticos N. ,N,,M- e M, isto é, os que são obti- 


dos directamente dos deslocamentos. A partir do 
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QUADRO I 


Pontos e 6 % Ka 5 Kp0 Forças fictícias (t) 
TES e HER 

3 " (m) (ºC) (ºC) f; m 

o | 2 | 695| —178| —603 | —32216145 | —63381308 0 11,0 

ol «4 |y7v50|—-2z42| —455 | —47265625 | —55847168 o 17,5 

o | 6 | 6,53|-—-491| —263 | —83495571 | — 24170579 0 12,0 

o| s | 384 | —-8,89 | -0,66 | -—88900000 | — 2088960 O |-258 

2 2 | 697 | —-212| —5,61| — 38480209 — 58786263 | 6,1 9,8 

2/4 | 744) —261| —438 | —50568750 | —52613885 | 5,6 17,1 

2 6 | 6,49 — 4,98 | —2,58 | — 84167189 | — 23729119 | 1,6 11,2 

21 8 |! 384 | —-8,90 | —0,65 | —88536461 |! — 2056490 |—-0,7 |-—26,4 

4 2 | 6,92 | —2,70 | —4,90 | — 48656249 — 51058042 | 8,9 11,0 

à | a 742 | = 2,25 3,97 | — 59704168 | —43864738 | 8,6 16,6 

4 | 6 | 645 |—513 | —2,44 | —86167969 ! —22128413 | 3,4 9,3 

al 8 | 379|—-8,91! —0,63 | —87939847 | — 2008418 —-11 | —29,0 

o 2 | 673|—3,35 | —4,29 | —58712241 | —42336289 | 10,0 | —124] 

6 | 4 | 678 | —3,91 — 3,52 | —69035938 | — 35170614 | 9,3 16,5 —21,7 

6. 6 | 596 | —582 | —2,24 | —90331248 | —17458620 | 3,8 

| 6 | 8 | 370 | —-901 | —0,58 | —86815107 | — 1758785 |-2,2 |—348 —18,3 

8 2 | 647 | —-408 | —3,69 | —68743749 | —33502347 | 9,4 

8 | 4 | 636 | —-47 | —3,05 | —78009375 | —26797860 | 84 15,0 | — 22,5 

8 | 6 | 5,56 | —6,43 | —1,98 | —93101044 | —13428095 | 20 2,8| — 24,5 

sl s | 3,55| —912| —0,51 | —84312503 | —1 359254 |—-2,9 | —42,9| — 18,2 
10 | 2 | 609 | —-5,05 | —2,99 | —78906250 | —22968750 | 8,8 12:86) 29 
10 | 4 | 5772 |—-5,79 | —2,47 | —87000783 | —17926001 | 6,2 419] —233 
10 | 6 | 501| —7,24 | —1.61 | —94459375 | — 8732743 |-0,3 |— 3,2] — 26,2 
10 | 8 | 3,34 | —9,32| —0o,41 | —81064586 | — 970305 |-37 | —56,2|] —19,7 
12 | 2 | 548 | —618 | —2,36 | —88193748 | —15356691 | 5,3 10,3| — 13,8 
12 | 4 | 517 | —-6,81| —1,94 | —91686721 | —11395746 | 2,0 41] 21,2 
12 | 6 | 4,46 | —7,96 | —1,26 | —92452081 | — 5390762 | 0,3 |—10,4| —28,1 
12 | 8 | 3,09|]—-9,51| —0o,28 | —76525781 | — 576243 |-3,3 |—72,1| —20,5 


valor dos primeiros calcularam-se as tensões 
7. e 3 ao longo do perfil central (fig. 10). 


Pela análise do andamento destas curvas veri- 
fica-se que para o diagrama de variações tér- 
micas considerado no exemplo se geram ten- 
sões de tracção a meio do paramento de jusante 
que atingem valores da ordem dos 20 kg' cm”. 
Estas tensões de tracção horizontal na secção de 
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fecho de jusante invertem-se próximo do coroa- 
mento atingindo no bordo livre valores máximos 
de compressão da ordem dos 30 kg/ cmº. 

No paramento de montante as tensões hori- 
zontais na secção do fecho da abóbado são quase 
exclusivamente de compressão, 

As tensões de tracção verticais na secção estu- 
dada são máximas no pé de jusante, onde atin- 
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gem valores máximos da ordem dos 20 kg/ cm”. 
No pé de montante as tensões verticais são de 
compressão, mantendo-se com o mesmo sinal 
numa extensa zona da secção de fecho, mas inver- 
tendo-se próximo do coroamento. 

O bordo livre causa uma importante pertur- 
bação no andamento das tensões horizontais e 
verticais reinantes nos paramentos de montante 
e jusante da secção do fecho. É de notar 


Cs um 7 dee vaca 


a importância das tensões resultantes das varia- 
ções térmicas que atingem, para o diagrama 
de oscilações térmicas analisado no exemplo, 
tensões de tracção e compressão da ordem de 1 3 
dos valores das tensões máximas de compressão 
na abóbada devidos à pressão hidrostática. 
Embora no exemplo estudado as tensões de 
tracção máximas se verifiquem em zonas onde 
reinam tensões de compressão devida à pressão 
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Material fabricado nas nossas oficinas, pronto para expedição 


Fabricação de aparelhagem eléctrica de alta e baixa tensão 
Equipamentos electromecânicos completos para centrais e fábricas 
Postos de transformação tipo monobloco 
Quadros blindados, capsulados e em armário 
Armaduras de iluminação de todos os tipos 
Contadores eléctricos e instrumentos de medida 
Telemedida e telecomando 

Instalações elevatórias de água 

Instalações de saneamento e esgotos 
Ferramentas de corte 

Motores Diesel industriais e ii 
Fogões eléctricos 

Frigoríficos 

Ventilação e condicionamento de ar 
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hidrostática, seria possível escolher diagramas de 
variações de temperatura que originassem ten- 
sões de tracção que se sobrepusessem às da 
pressão hidrostática. Só um estudo analítico 
e sistemático permitiria esclarecer devidamente 
a importância das tensões térmicas e a sua dis- 
tribuição na abóbada para diversos diagra- 
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mas-tipo. Este problema tem sido abordado 
por via experimental [5], [6], mas aí as 
técnicas de ensaio sobre modelo complicam-se 
de tal modo que se acha preferível numa 
base preliminar, para apoio e esclarecimento 
dos resultados experimentais, o seu estudo ana- 
lítico. 
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DISJUNTORES 
TRIPOLARES 
R 920 


Equipados sempre com reles 
bimetálicos e daisparadores 
electromagneticos, são os apa- 
relhos indicados para à pro- 
tecção de qualquer equipa- 
mento eléctrico trifásico 


ARRANCADORES 
ESTRELA - TRIÂNGULO 
K 987 


Dotados com contaciores de 
ruplura no ar, permitem a liga- 
ção de motores trifásicos de 
rotor em curto-circuito até 210 
Kw, 500 V, 50 €;s 


Peça catálogos e informações à 


SIEMENS 


Para 
proteger 
e 
comandar 
motores 


CONTACTORES 
K 915 


Caracterizam-se essencialmen- 
te pela elevada frequência de 
manobras, alto poder de corte 
e longa duração, podendo ser 
fornecidos com relés bimetá- 
licos 


INVERSORES 
AUTOMÁTICOS 
DE MARCHA 


Utilizando os contactores K 915, 
são próprios para uma potên- 
cia até 35 KV em 500 V, poden- 
do tambem ser equipados com 
relés térmicos, É 


lv 17 66 096 


vale 

a pena 
consultar a 
Siemens 


SIEMENS - COMPANHIA DE ELECTRICIDADE, S.A.R.L. 


Sede: LISBOA -1 — Avenida Almirante Reis, 65 — Tel 536921/0 
Delegação Técnica: PORTO — Rua das Carmelitas, 12 — Tel 28943/7 
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CORRENTES DE 
TRANSMISSÃO 


De :” até 5” segundo as normas 


- (o BS, DIN, AFNOR, ASA e API. 
dio (É 
am WI RODAS EM DEPOSITO 
tt) para uso com correntes RENOLD 
eU) ni de transmissão e de transporta- 


produtos : 


RENOLD ami 


organização 
mundial no fabrico 
de correntes 

de precisão 


a 
E 
a SA , 


CORRENTES NORMAIS PARA 
TRANSPORTADOR 

com adaptações para todos os fins. 
Gama 1.360 kg (3.000 libras) a 
38.560 kg (85.000 libras). 


UNIOES DE ACOPLAMENTO 
por Corrente, Cruzeta de Borracha 
e Disco Flexivel até 2.500 CV. 


— 


EMBRAIAGENS 
para rotação livre, movimento 
divisor e sem retrocesso. 


K 


ll, 
II 


Para informações ou fornecimentos, escreva, visite ou telefone: 


HARKER, SUMNER & Ga., Lda. 


PORTO: 38, Rua de Ceuta, 48 
Apartado No. 162 
Tel: 27054-P.P.C. (4 linhas) 


LISBOA-2: 14, Largo do Corpo Santo, 18 
Tel: 324823-35124 
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O DIMENSIONAMENTO E A OPTIMIZAÇÃO DE PAVIMENTOS 
PRÉ-FABRICADOS (*) 


RESUMO 


Definem-se neste trabalho os esquemas de dimensiona- 
mento analítico dos pavimentos pré-fabricados com vigotas 
de betão pré-esforçado, utilizados no LNEC, e discutem-se 
os parâmetros que estão na base da optimização das vigotas, 
na fase do seu dimensicnamento, 

Finalmente, analisam-se diversos critérios de optimi- 
zação duma gama de pavimentos daquele tipc, com vista à 
sua redução aos elementos de aplicação mais favorável, 


1 — INTRODUÇÃO 


Uma das formas de pré-fabricação ligeira em 
que mais se tem apoiado a construção de edifi- 
cios do tipo tradicional tem sido a utilização de 
pavimentos pré-fabricados. Em Portugal verifi- 
cou-se, na última década, um notável incremento 
na aplicação desses pavimentos constituídos com 
vigotas de betão ou elementos cerâmicos, pré- 
-esforçados, assistindo-se a uma multiplicação 
do número das respectivas unidades de fabrico. 

Por força da legislação portuguesa, relativa às 
edificações urbanas, os novos materiais e pro- 
cessos de construção a aplicar no país, para os 
quais não exista uma experiência significativa de 
uso, devem ser apreciados e homologados pelo 
Laboratório Nacional de Engenharia Civil (LNEC). 
Como consequência directa deste facto, os pavi- 
mentos não tradicionais de betão armado ou pré- 
-esforçado — ou ainda semelhantes, do ponto de 
vista de dimensionamento e funcionamento, como 
os de cerâmica armada ou pré-esforçada — têm 
sido submetidos à homologação do LNEC, sendo 
incluída, nos planos de apreciação, a realização 
de cálculos de verificação de dimensionamento, 
efectuados por via automática, no computador 
NCR Elliot 803-B, de acordo com as dispo- 
sições da Circular da Informação Técnica (C..T.) 


por J. VASCONCELOS PAIVA 


Engenheiro Estagiário para Especialista 
Serviço de Edifícios e Pontes 
LNEC — Lisboa 


SYNOPSIS 


The analytic design achemes in use at the Laboratório 
Nucional de Engenharia Civil for prefabricated floora 
with prestressed concrete joists are presented, and the para- 
meters which govern the optimization of joists im the design 
stage are discussed. 

Finally several optimization criteria in use for diffe- 
rent flocrs of this type are analysed qith a view to selee- 
tiny the floors with optimum application. 


n.º 29 — “Pavimentos constituídos por vigotas de 
betão pré-esforçado. Regras para o seu dimensio- 
namento analítico e execução” [1]. 

Graças à experiência colhida no LNEC com a 
realização dos cálculos de vários milhares de 
pavimentos deste tipo, até esta data, foi possível 
estabelecer algumas vias de dimensionamento 
das vigotas, e definir os critérios que conduzem 
à optimização destas, e, paralelamente, das 
gamas de pavimentos que as integram. 

O objectivo deste trabalho consiste precisa- 
mente na definição dessas vias e dos critérios de 
optimização a seguir, quer para as vigotas, na 
fase do seu dimensionamento, quer para a gama 
de pavimentos constituída com base naquelas, a 
qual se procura reduzir aos elementos de utili- 
zação mais favorável. 


2 — DIMENSIONAMENTO ANALÍTICO DOS 
PAVIMENTOS PRÉ-FABRICADOS 


2.1 — Constituição duma gama de pavimentos 


Uma gama de pavimentos pré-fabricados com 
vigotas de betão pré-esforçado constitui-se com 
uma série de vigotas, tendo em geral a mesma 
forma de secção de betão e diferentes armaduras 


(*) Comunicação apresentada às II Jornadas Luso-Brasileiras de Engenharia Civil — Rio de Janeiro e S. Paulo, 1967 
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de pré-esforço, e um conjunto de blocos de 
cofragem, de diversas dimensões. 

Da associação em obra de cada bloco de 
cotragem com os vários tipos de vigotas (desde 
que se verifique em todas a igualdade de secções 
de betão), e podendo ser variável a altura da 
lâmina de compressão de betão complementar, 
que solidariza entre si aqueles elementos pré- 
-fabricados, resultam diferentes pavimentos per- 
tencentes ao mesmo tipo. 

Um tipo de pavimento é, por conseguinte, 
definido pela utilização de um bloco e de vigotas 
com forma idêntica, sendo a distinção entre 
pavimentos do mesmo tipo conferida pela arma- 
dura da vigota utilizada, e pela altura da lâmina 
de compressão, normalmente de 3, 4, 5 ou 6 cm. 

Os tipos de pavimentos podem ainda ser 
simples ou duplos, consoante as nervuras resis- 
tentes entre blocos são constituídas por uma 
vigota, ou por duas vigotas colocadas a par. 

No caso geral, e a menos que haja que prever 
soluções alternativas — por exemplo, quando cer- 
tas carências no mercado dos aços aconselhem a 
utilização de vigotas sensivelmente equivalentes 
do ponto de vista resistente, mas armadas com 
varões de diferentes diâmetros — bastarão quatro 
ou cinco tipos de vigotas, se convenientemente 
estudados, para satisfazerem as necessidades 
duma gama de pavimentos, entre limites bas- 
tante largos. 

É também importante a escolha criteriosa das 
dimensões exteriores dos blocos de cofragem a 
utilizar, das quais vão depender directamente o 


afastamento entre nervuras e a altura dos diversos 
pavimentos, com vista a obter-se um elevado 
rendimento e um bom escalonamento da gama. 


2.2 — Determinação da capacidade resistente 
dos pavimentos 


Além do momento de inércia da secção trans- 
versal dum pavimento, e dos módulos de elastici- 
dade dos materiais que o constituem — grandezas 
responsáveis pela deformabilidade do pavimento, 
que não deve ultrapassar certos valores limites [1] 
— a capacidade resistente desse pavimento fica 
definida pelos valores do seu momento flector de 
serviço, M., e do seu esforço transverso de ser- 
viço, T,. O pavimento poderá suportar determi- 
nadas condições de vão e de solicitação (do tipo 
I), se os esforços actuantes M e T, originados 
nessas condições, obedecerem às seguintes rela- 
ções evidentes: 

M<M, 
TT 


Vejamos então como se calculam os esforços 
de serviço, para analisarmos em seguida como é 
possível actuar no dimensionamento, para me- 
lhorar os valores desses esforços, em especial o 
momento M,. 

Antes disso, torna-se necessário que se 
indiquem quais as tensões permitidas nos mate- 
riais, de acordo com a C.I.T. n.º 29. 

Conforme a classe de betão das vigotas, B350 
ou B450, as tensões admissíveis nestes elementos 
em fase elástica são dadas no quadro seguinte: 


Cr E TS 


Tipos de esforços 


Compressão resultante da flexão  — - 


Corte 


Tensões admissíveis 


kef'cm? 

== 350 | B 450 
mid [HE] —m% 
tenso es — 110 — 140 
na fibra inferior u + 10 
na fibra superior = +30 o 
sem armaduras sua esforços ines E E 

versos | 
“com armaduras para mini ra o » 


versos 
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